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ABSTimCT 


The  general  method  of  Oorodniteyn  (14),  for  integrating  first-order  partial 
differential  equations  of  mixed  type,  is  applied  to  calculate  the  flow  of  a  uniform 
supersonic  or  hypersonic  stream  past  a  circular  cylinder.  Using  polar  coordi¬ 
nates  referred  to  the  center  of  the  cylinder,  the  equations  of  motions  are  reduced 
to  an  approximating  system  of  ordinary  differential  equations  in  the  angular 
coordinate,  by  prescribing  the  variation  of  key  dependent  variables  with  radial 
distance.  In  the  first  approximation  this  variation  is  linear,  in  the  second 
approximation  it  is  quadratic  and  so  on.  The  boundary  conditions  applied  to  the 
ordinary  differential  equations  are  given  by  symmetry  conditions  and  conditions 
that  certain  equations  pass  smoothly  through  singularities  near  the  sonic  line. 

Flow  of  a  perfect  gas  past  a  circular  cylinder  is  calculated  for  free-stream 
Mach  numbers  of  2. 13,  2. 5,  and  4.  0  using  the  first  and  second  approximations, 
end  for  free-stream  Mach  numbers  of  3.  0  and  5.  0  using  the  first,  second,  and 
third  approximations. 

Results  show  excellent  agreement  with  experimental  data.  The  method  can  be 
extended  to  calculate  the  flow  of  a  nonperfect  gas  past  cylinders  or  bodies  of 
revolution  of  general  shape. 
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1  -  shock  wave 

2  -  circular  cylinder 

2  Subdivision  of  the  region  of  the  integration  in  the 
various  approximations: 
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3  Shock  wavei  sonic  line,  and  the  characteristics  that  bound  the 
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5  Pressure  distribution  on  the  surface  of  a  circular  cylinder 
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7  Pressure  on  the  surface  of  the  circular  cylinder  in 
various  approximations  at  >  3-0 

— N  -  1; - -  N  *  2;  AAA  N  -  3 
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14  Comparison  of  the  distribution  of  pressure  along 
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■ - -  ■  N  ■  2  ;  #  #  #  experiment 
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point) 

Mach  number 
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Y 
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k 
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a 
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SYMBOLS  USED  (Cont'd) 


H  >  kp  ■«-  pu^ 
h  ■  rv  . 

Subscripts: 

■*  values  ahead  of  the  shock  wave 

i  values  on  the  shock  wave 

i  ••  2,  3. . .  .  .  N  values  on  the  intermediate  lines  behind  the  shock  wave 

0  values  on  the  surface  of  the  body. 


INTRODUCTION 


Although  the  problem  of  flow  past  a  body  with  a  detached  shock  wave  has  been 
investigated  by  many  authors,  most  of  them  only  considered  an  approximate 
statement  of  the  probleni<  The  difficulty  in  problems  of  this  type  is  that  the 
gas  flow  is  of  mixed  type,  where  on  one  side  of  the  transition  (sonic)  line  the 
motion  is  subsonic  (elliptic  system  of  equations),  while  on  the  other  side  the 
speed  is  supersonic  (hyperbolic  system  of  equations).  The  velocity  varies 
over  a  rather  wide  range,  from  zero  (at  the  stagnation  point)  to  supersonic, 
and  ordinary  linearization  methods  fail.  In  addition,  the  boundary  conditions 
of  the  problem  are  specified  both  on  the  body  and  on  the  shock  wave,  but  the 
shape  and  position  of  the  latter  are  not  known  is  advance. 

In  many  papers,  the  solution  of  this  problem  is  given  in  series  form.  Thus, 
G.S.  Smurov  (1)  in  an  analysis  of  the  flow  of  a  supersonic  stream  of  gas  past 
a  circular  cylinder  writes  the  stream  function  in  the  vicinity  of  the  symmetry 
axis  in  the  form  of  a  series,  coefficients  of  which  are  determined  from  the 
boundary  conditions  of  the  problem.  Smurov  finds  the  distance  ^  from  the  sur¬ 
face  of  the  body  to  the  shock  wave  along  the  axis  of  symmetry,  and  also  the 
pressure  distribution  on  the  body  near  the  stagnation  point.  His  calculated 
data  are  presented  in  Fig.  11,  and  it  is  seen  that  they  differ  considerably  from 
experimental  results.  It  must  be  noted  that  Smurov* s  paper  is  one  of  the  first 
in  this  field.  I.  A.  Panichkin  (2)  and  H.  Melkus  (3)  expand  the  velocity  along 
the  axis  of  symmetry  in  a  ieries  and,  assuming  the  shock  wave  to  be  in  the 

form  of  a  hyperbola,  find  the  ratio - for  various  values  of  the  incident  stream 

*0 

Mach  number  M„(fg  is  the  radius  of  curvature  of  the  wave  at  its  crest).  C.C. 
Lin  and  S.I.  Rubinov  (4),  also  using  a  series  es^ansion,  obtain  a  relation  be¬ 
tween  the  position  of  the  shock  wave  and  the  radius  of  curvature  of  the  body 
at  the  nose.  The  method  of  series  expansion  makes  it  possible  to  investigate 
only  a  small  region  of  the  flow  in  the  vicinity  of  the  axis  of  symmetry  or  else 
in  the  vicinity  of  the  stagnation  point  of  the  body,  and  since  there  is  no  inform¬ 
ation  on  the  region  of  convergence  of  these  series,  attempts  by  several  authors 
to  extend  these  solutions  to  the  sonic  line  must  be  deemed  unsuccessful. 

The  possibility  of  obtaining  an  nalytic  solution  behind  a  curved  shock  wave  in 
the  hodograph  plane  has  also  bi-.i)n  investigated.  Q.  Guderley  (5)  and  A. 
Busemann  (6)  considered  the  ci^se  vhen  the  Mach  number  of  the  incident  stream 
is  close  to  unity.  This  enabled  them  to  assume  the  flow  behind  the  wave  to  be 
isentropic,  and  so  to  use  the  Chaplygin  differential  equations.  K.  Hida  (7), 
using  two  particular  solutions  of  Tricomi's  equation,  obtains  formulas  to  esti¬ 
mate  the  asymptotic  position  of  the  detached  shock  wave  for  flow  of  a  gas  with 
M„  close  to  unity  past  blunt-nosed  bodies. 
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The  flow  in  the  vicinity  of  the  crest  of  the  wave  is  consiuered  by  Hida  to  be 
i»entropic>  He  considers  this  flow  in  the  hodograph  plane  and  neglects  the  in¬ 
fluence  of  the  streamline  corresponding  to  the  surface  of  the  body,  in  view  of 
the  closeness  of  to  unity.  However  in  the  general  case,  one  cannot  neglect  the 
influence  of  the  vorticity  behind  the  detached  shock  wave.  It  is  difficult  to  ob-  . 
tain  an  analytic  solution  for  rotational  flow  in  the  hodograph  plane  because  it 
is  difficult  to  write  down  the  boundary  conditions  along  the.  surface  of  the  body, 
and  also  to  transform  these  streams  into  the  hodograph  plane. 

K.  Tamada  (8)  and  K.  Hida  (9)  found  the  position  and  shape  of  the  shock  wave 
for  flow  around  a  sphere  and  a  circular  cylinder  by  assuming  that  the  wave  is 
followed  by  rotational  flow  of  an  incompressible  liquid.  T.  Kawamura  (10)  ex¬ 
pands  in  this  case  all  the  hydrodynamic  variables  in  a  series  in  the  vicinity  of 
the  stagnation  point  and,  assuming  a  hyperbolic  or  a  parabolic  shape  of  shock 
wave,  he  determines  th^  distance  from  the  wave  to  the  surface  of  the  cylinder 
or  sphere.  The  results  obtained  by  these  authors  are  in  fair  agreement  v/ith 
experiment  for  flow  around  a  sphere,  but  the  agreement  is  much  worse  for  flow 
around  a  circular  cylinder,  as  can  be  seen  from  Fig.  11.  S.  Uchida  and  M- 
Yasuhara  (11)  propose  a)',  approximate  graphical  me.  >od  for  calculating  the 
flow  behind  a  curved  shock  wave,  in  which  they  must  first  specify  the  com¬ 
plete  flow  pattern.  The  .boundary  conditions  on  the  wave  are  satisfied  approxi¬ 
mately,  and  the  accuracy  of  the  method  diminishes  rapidly  on  approaching  the 
sonic  line.  The  distance  cq  from  the  body  to  the  wave,  obtained  by  the  authors 
at  If.  >  2  0,  is  in  good  agreement  with  experiment. 

Another  way  of  obtaining  an  approximate  solution  as  proposed  by  J.  W.  Maccoll 
and  J.  Codd  (12)  and  by  A.R.  Mitchell  (13);  specifying  the  experimentally  ob¬ 
tained  shape  and  position  of  a  shock  wave,  they  used  the  relaxation  method  to 
calculate  the  flow  around  a  body  with  a  flat  nose.  J.  W.  Maccoll  and  J.  Codd  con¬ 
sidered  the  flow  behind  the  wave  to  be  free  of  vorticity  (  m.  >  1.} ). 

The  use  of  electronic  computers  makes  it  possible  to  obtain  complete  results 
of  any  required  accuracy  for  the  problem  In  its  rigorous  form.  A  sufficiently 
effective  method  for  use  v^th  electronic  computers  is  one  in  which  integration 
of  the  syi^tom  of  nonlinear difference  equations  is  reduced  to  a  numerical 
solution  of  a  certain  i^ptroximating  system  oi  ordinary  differential  equations. 
This  method  was  proposed  by  Academician  A.  A.  Dorodnituyn,  and  was  reported 
together  with  certain  results  of  the  calculations  at  the  Third  All-Union  Mathe¬ 
matical  Congress  (14). 

In  another  paper  (15),  We  ua«ed  this  method  to  obtain  a  solution  of  the  problem 
of  the  flow  past  a  plane  body  of  arbitrary  shape  with  an  axis  of  symmetry  and  a 
detached  shock  wave.  In  the  present  paper,  we  give  the  formulas  and  results  for 
a  circular  cylinder. 
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All  the  calculations  were  made  on  the  BESM  high-speed  eiuctronic  computer 
of  the  USSR  Academy  of  Sciences,  built  at  the  Institute  for  Precision  Mechanics 
and  Computational  Engineering  under  the  leadership  of  Academician  S>  A*  Lebedev. 

I  express  my  deep  gratitude  to  Comrades  E.S*  Bogomolova,  N-S>  Malinina, 

V.  P.  Blasova,  N<V.  Mel'tsis,  A- !•  Bykova,  L.  V.  Papandina,  A.N.  Belova, 
and  others  who  helped  in  checking  the  formulas,  in  the  performance  of  the  pre¬ 
liminary  calculations,  and  also  in  the  processing  of  the  data  obtained. 


I.  STATEMENT  OF  THE  PROBLEM 


Let  a  plane -parallel  supersonic  stream  1 )  of  an  ideal  gas  be  incident  with 
constant  velocity  on  a  circular  cylinder.  A  shock  wave,  the  shape  and  posi¬ 
tion  of  which  -are  unknown  beforehand,  is  produced  in  front  of  the  cylinder.  It 
is  required  to  calculate  the  mixed  rotational  flow  of  the  compressed  gas  in  the 
region  bounded  by  the  shock  wave,  the  axis  of  symmetry,  the  surface  of  the 
body,  and  the  first  limiting  characteristic  (or  characteristics)  passing  between 
the  wave  and  the  body  and  bounding  the  minimum  region  of  influence.  Outside 
this  region,  the  calculation  of  the  flow  can  be  made,  say,  by  the  method  of 
characteristics. 

We  introduce  dimensionless  quantities  which  relate  the  velocity  to  the  maximum 
velocity,  the  pressure  and  density  to  the  stagnation  pressure  and  density  ahead 
of  the  shock  wave,  and  linear  dimensions  to  the  radius  of  the  cylinder.  The 
equations  of  motion,  continuity,  and  energy  (adiabatic  condition)  then  become, 
in  vector  form 


cutl'^  X  V  + 


v(fcp) 

P 


0 


V(^)  -  0 


PY 


(1.1) 


where  w,  p  ,  and  p  are  the  velocity,  pressure,  and  density  behind  the  shock 
wave,  and  y  is  the  adiabatic  exponent. 


(for  air  v  ■  1,40  andk  ■  1/V). 

The  dimensionless  form  of  the  Bernoulli  integral  is 

-  1  .  (1.2) 

P 

The  equation  of  continuity  leads  to  the  conclusion  that  there  exists  a  stream 
function  such  that 


di/t 

«k 


dt/i 

pr.-.-rpu 


(1.3) 
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hence 


(1.4) 


dr 

If  we  introduce  a  vorticity  function  ^ ••p/^Xfentropy  s  =  C^ln^),  we  obtain  from 
the  adiabatic  condition  and  from  (1.  3) 

^  S”  -  ^(Vlr)  (1.5) 


after  which,  taking  (1.2)  into  account,  we  find 

1 


P  -  P 


(1  -  w^)  -  ^  y  ~  1 


p  - 


where 


r-(l«w2)y“l  . 

The  velocity  of  sound  is  defined  as  follows: 
■)  kp  y  -  1 

-  y  — -  - -  (1  -w2)  . 


We  rewrite  the  system  (1.  1)  in  polar  coordinates  r  •  9  (Fig.  1),  and  introduce 
the  Bernoulli  integral  in  place  of  one  of  the  equations  of  motion.  We  then 
obtain  a  system  equivalent  to  (1. 1).  Together  with  Eq.  (1.4),  the  new  system 
will  be 


diH  d% 

IT  *  99 


dtH  di 
dt  *  99 


_-p(r__«,  .  I  (1.10) 

Here  H  -  kp  pu^,  s  -  pur,  g  »  kp  +  p  h  -  ri»,  t  -  rr,  and  p  ,  p  ,  r  are  found 
from  (1.6)  --  (1.8).  The  unknown  function's  in  this  system  are  u  ,  t  ,  ^,  and 
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We  now  write  down  the  boundary  conditions  of  the  problem.  On  the  surface  of 
the  cylinde  r  (  r  =  rg  =  1  )  we  have 


/2ziy_L 

y2_,  \y+l/  4y 

We  write  the  equation  of  the  wave  in  the  following  form 


1(0)  -!+<(«)  A1 

where  e(0)  is  the  distance  from  the  surface  of  the  body  to  the  wave  along  the 
ray  0  =  const.  The  relations  for  the  tangential  and  normal  components  of  the 
velocity  on  the  shock  wave  are,  according  to  Sauer  (16), 


Furthermore 


*  .  y 

- and  w_  ■  -  - 


where  a  is  the  angle  of  inclination  of  the  wave  to  the  direction  of  the  incident 
stream,  and  the  index  «  denotes  values  ahead  of  the  shock  wave. 

Hence 

•-  .  /  1  \  1 
w_  ■  — —  MS  2e  I  1  -  — '  '  —  ■  1  / 

'  V  I 

w_  »  1  -  — —  sin^  oil  —  - i -  y  I  \  (1.12 

‘  \  m2  )y  I 

The  necessary  conditions  for  the  veldcities  are  then  written 


u  >■  Wy  sin  0  —  cos  0  i 
T  »  Wjj  sin  +  w„  cos  0  .  J 
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II.  METHOD  OF  SOLUTION 


As  already  noted  iii  the  introduction,  the  problem  posed  is  solved  by  the  method 
of  integral  relations  (14). 


We  set  up  an  approximating  system  for  the  system  ot  equations  (1.  10).  For 
this  purpose,  w'e  draw  between  the  body  and  the  wave  a  set  of  N-1  lines  that  are 
equidistant  along  r  : 


r  -  rj  ( «) 


N-i  \ 

N  ) 


B1 


These  lines  break  up  the  region  of  integration  into  N  strips  (see  Fig.  1).  We 
denote  all  the  quantities  on  the  intermediate  line  by  the  index  i  ,  on  the 
shock  wave  (i  =  1 )  by  the  index  1,  and  on  the  surface  of  the  body  by  the  index  0. 


By  integrating  the  first  two  equations  of  (1. 10)  along  an  arbitrary  ray  6  =  const 
from  the  surface  of  the  body,  t  s  1  ,  to  the  boundary  of  each  of  the  strips,  we 
obtain  2N  independent  integral  equations 


J  fjHi  -  Il0  •  ^  g(r,  d)dr 


dr-., ho-  0 
I 

(i  -  1.2 . M)  .  (2.2) 

We  now  approximate  any  integrand  f(r,  d  )  by  means  of  an  interpolation  poly> 
nomial  of  degree  N  in  r  ,  taking  the  interpolation  points  to  be  the  boundaries  of 
the  strips 


(2.3) 


where  the  coefficients  (d)  depend  linearly  on  the  values  of  the  corresponding 

functions  on  the  boundaries  of  the  strips. 

« 

Next,  writing  down  the  two  last  equations  of  (1. 10)  along  each  of  the  N-1  inter¬ 
mediate  lines,  r  -  r,  (d)  (  tft  and  4  on  the  wave  and  on  the  body  are  determined 
from  the  boundary  conditions),  and  taking  Eq.  (1. 16)  into  account,  we  obtain 
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an  approximating  system  consisting  of  4N  - 1  equations  for  the  unknowns  c  ,  a  , 
tq  •  uj  ,  vj  ,  I  ( i  2,  3.  •  •  ■  .  N ).  In  normal  form,  this  system  can  be 
written  schematically  as: 


—  -  -  (1  +  «)C0t(ff  +  ®) 
dv 


dd  y  -  1 
y+  1 


d^i  Pdfj  -j 

du;  I  — -  s:  d  In  ^  d^: 

t.  6:  "  ^  ^ 


y  -  1  +  2u? 


y  +  1 


(^i) 


i  -  2.  3,  ....  N  . 


1-  w2 


k  [•■  -  '‘(i  -  “?)] 
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»?)  J  Gl 


"  “7"  -  "i 

_  1 


In  the  tyatem  (2.4),  f  •  Eq  •  £>  •  m(*  andq'  are  holomorphic  functions  of  6  and 
of  the  quantitiea  required,  ano  are  defined  in  the  range  of  integration  (the  form 
of  theee  functions  depends  on  N  ),  and 


**  **  ^  ff. 

4^1  49  40  d^l  ^ 

where  die  de  is  determined  from  (1. 14),  d^ij/  40  is  determined  from 
(1.15),  and  49  /  40  *  f  is  determined  Irom  (2.4). 


We  have  thus  obtained  a  system  consisting  of  3N  ordinary  differential  equations 
and N-1  finite  relations.  All  the  boundary  conditions  on  the  body  and  on  the 
shock  wave  are  automatically  satisfied,  as  can  be  seen  from  the  scheme  by 
which  the  approximating  system  has  been  constructed. 


The  integration  of  the  resulting  system  is  carried  out  numerically  beginning  at 
the  axis  of  symmetry^  -  0  ,  where  <r  -  s/2,  tq  -  rj  «  -  0  ,  while  the  N 
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initial  values  of  the  fvmctions  <q  and  U|  ( i  >  2,  3, _ _  M  )  are  unknown  (we  have 

N  parameters  at  9  -  0  ).  The  problem  of  supersonic  flow  past  bodies  with  a 
detached  shock  wave  is  of  mixed  type.  As  is  well  known,  the  distrubances  that 
arise  in  the  supersonic  region  past  the  boundary  characteristic  are  not  propagated 
in  such  problems  into  the  subsonic  region.  We  can  therefore  obtain  no  additional 
conditions  for  the  determination  of  the  remaining  initial  values  of  the  functions 
required  on  the  upper  boundary  of  the  region  of  integration.  This  fact  is  specifi¬ 
cally  reflected  in  the  approximating  system.  It  is  seen  from  its  structure  that  in 

the  vicinity  of  the  sonic  line  |  ■>  — - )the  N  equations  of  the  system  will  have 

V  y  +  1  / 

floating  singular  points,  and  special  conditions  must  be  satisfied  if  a  continuous 
transition  through  these  points  is  to  be  possible. 


Actually,  at  points  where  w~ 


_  ^  ~  ^  and  wf  - 
0  y  +  1 

the  set  of  these  points  the  "singular  line"),  we 


y  —  1  2u;^ 

-  +  — —  (we  shall  call 

y  +  1  y+1 

should  have  respectively 


Eg  -  0  and  Ej  -  0  (i  -  2,  3 . N)  .  (2.5) 

Otherwise  the  derivatives 

<l»0  ‘•''2  ‘*’^3 

dT  ’  dT  ’  dT . IT 

become  infinite  on  the  singular  line;  i.  e, ,  the  accelerations  become  infinite, 
and  the  motion  cannot  continue  past  this  line  —  the  Singular  line  itself  will  be 
the  limiting  line,  and  all  solutions  will  lose  their  physical  meaning. 

To  satisfy  the  N  Conditions  (2.  5),  we  have  at  our  disposal  N  parameters  at 
0  m  0  ,  Thus,  the  requirement  that  the  motion  be  continuous  on  the  singular 
line  indeed  supplies  us  the  conditions  necessary  to  make  the  problem  definite. 

In  reference  15,  we  investigated  the  character  of  the  behavior  oi!  integral  curves 
in  the  vicinity  of  the  singular  points  (the  singularities  were  of  the  "saddle"  type). 
We  also  showed  there  that  the  approximating  system  (2,4)  has  a  unique  solution, 
is  holomorphic  in  the  entire  region  of  integration,  and  satisfies  the  conditions 
at  9  -  Oand  on  the  singular  line. 

We  present  computational  formulas  for  N  -  1,2, 3. 

First  approximation  (  N  "  1 ).  There  are  no  intermediate  lines  between  the  body 
and  the  wave  (Fig.  2ci). 
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The  approximation  is  of  the  form 


-  to  B6 

where  (9)  and  (ff)  are  the  values  of  f  (r,  9)  on  the  cylinder  and  on  the 
wave,  respectively.  The  unknowns  are  the  functions  e  ,  a  ,  tq  ,  and  the 
functions  that  remain  to  be  determined  on  the  body  and  on  the  wave  are  found 
from  the  boundary  conuitions  (1. 11)  and  (1. 13  -  1.  15).  The  approximating 
system  is  readily  obtained  from  (2.4)  (first  three  equations).  In  this  system 


*1  "  ~  1^*1  ^  2  (  1  + 1)  flj J  +  *0  +  *1 

,  1  d*  do  /  2  \ 

*0  -  *0)7^  -  Gl  -jj-  +— jhl  • 


When  9  m  0  we  have  tq  ■  0  and  a  «  ir/2  ,  while  »(0)  «  ^  is  the  one  unknown 
parameter,  determined  from  the  following  condition:  we  should  have  Eg  -  0 

5  y  “  * 

at  the  point  Ag,  where  ^ " 

Second  approximation  ( N  ■  2  )«  We  introduce  one  intermediate  line  (Fig,  2b). 
The  approximations 


ttt.9)  -  ♦  (dfj  -  t^  -  3fo) 


>«.  -  •  <•>  [—]' 


where  f|(d)  andfj  (9)  are  the  values  off  (r,  9)  on  the  cylinder,  on 

the  wave,  and  or,  whe  central  line,  respectively.  The  unknown  functions  are 

<  .  f  "2  •  ’'2  *  ^2  •  ^  • 

In  the  system  (2. 4)  we  have 


-  (2  +  *)  H2J  I  ♦  ,1  — 


♦  ( 1  ♦  «)  H, 


»1  -  *0 
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•2  "IT  [«!  ^*  5«0  -  H,  -2{2+0H2]  +  y  +  Ij 

“o'-  [d+Ohi  -  (2  ♦  Ohjj  -  (to  +  3ti  -  4t2)^| 


da 

ti  -  G,  —  -  hj 


*1  1  r  dt 


—  +  (2  +  «)h2  -y(l  ♦  «)hi 


When  0  >  0  we  have  to  «  ®  •  *  “  1^/2  ,  and  ^  and  uj  (0)  are  the  two 

unknown  parameters,  which  are  determined  from  the  following  conditions:  at 
the  point  Aq  ,  where  t2  ,  (y  _  ^  should  have  Eq  -0  .  and  at  the  point 

Aj  ,  where  w2  -  (  y  _°i  +  2u2 )/(  y  +  1 ),  we  should  have  E2  -  0  . 

Third  approximation  (  N  ■a).  We  introduce  two  intermediate  lines  (Fig.  2c). 


The  approximations  are 


t(t.0)  -  to  +  j(2(,  -  U(o  -  9f2  ♦  »»<j) 

*7"'' [^fvj 

where  f.  (0  -  f  (tj.  9) .  i  -  0,  1,  2,  5- 
The  unknown  functions  are 

*2’ ^2’''3’ ’“3' ^2' ^3’ ^2' ^3  ‘ 

The  functions  appearing  in  the  right  halves  of  the  system  (2. 4)  are  of  the 


form 


•j  “  “  (24A2  —  13A|  —  3Aj) ,  *2  -  |J^^3  “  *’  ^^2^  ’  *3  *  ~  *^2  “  ^3^  * 

tQ  “  (24B2“13Bj-3B3-16t|'),tj  ■  Gj  — y  -  h|  ,  t^  ■  y  (0*2  “  ♦  *3  “  i0*i ) 
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tj  «  —  (7Bi  16 B2  +  Bj  +  8  )  ; 


Bl  -- 


3  (*2  *3^  “  ^  *1 


]  *  [<*  +  -  ”o]|-  <0  -  *1  -  5  C«2  +  83). 


-"Z  -  f  I  ^  ♦  t)  -  "»]  I  -  *»  -  '2  -  ■'»3  • 

*J  -f  |[•I  -  5  <*2  ♦  -J*]  -  H„  I  -  9,„  -  „  .  5,j 

*l  -  7  I  ["O  -  2>i  ♦  3  (tj  t  .,)]  ^  .  8  (1  *  ,)  k,  j 
»2  -  7  Ift  -  ♦  «3]  ^  t  hj  I  , 

»3-7{[»'t)**.  -»«I*',)]  1*32  (‘*f)'’3| 

When  6  -  0  we  have  vq  -  vj  "  ^3  -  lt'2  -  1^3  <■  0.  a  ■  n/1  »  *0  »  “2 

03(0)  are  the  three  unknown  parameters,  which  are  determined  from  the 
following  three  conditions: 


+  5*2  -  1983 


Rq  -  0  u  dw  poise  Aq,  wkois  -(.y  -  l)/(y  +  1) 

Ej]  -  0  M  Ike  poise  A21  »fc«eo  "^y-  1  +  2u^^/(y  +  1) 

E3  -  0,  se  (he  poise  A3,  ehoes  wj  -  ^y  -  1  ♦  20*^  /(y  +  1)  . 

To  conclude  this  section,  we  note  that  in  this  problem  it  was  found  most  con¬ 
venient  to  use  polar  coordinates  and  to  approximate  the  functions  in  terms  of 
r.  Also  we  have  attempted  to  approximate  the  functions  in  terms  of  0  ,  and  also 
introduced  the  coordinstee  4(4"  const  being  orthogonal  to  3^  «  const),  ^  ,  t  , 
..nd  ^  ,  iy(  ^  (  f  -  tQ  (6)  )  /  e(6),  17  -  6/  6j  (  e ),  where  ( r )  is  the  equa¬ 

tion  of  the  limiting  characteristic),  but  the  approximating  systems  obtained 
thereby  were  very  unwieldy  and  had  a  large  number  of  parameters.  Next,  if 
we  take  as  the  initial  system  of  the  differential  equations  a  system  equivalent 
to  (2. 4),  in  which  the  unknown  functions  are  u  ,  v  ,  p  ,  and  p  ,  we  obtain  an 
approximating  system  with  2N-1  parameters  at  9-0  and  2N-1  conditions  at  the 


singular  points.  During  the  formulation  of  the  problem  we  carried  out  the 
calculation  in  a  system  like  this,  also,  but  even  in  the  first  approximation 
( N  «=  1  )  the  results  of  the  calculation  were  quite  close  to  each  other  in  both 
systems.  For  example,  for  =5.0we  had  f  -  0.471  in  the  system  with  the 
unknown  functions  u  ,  v  ,  and  xp  ,  and(  •  0.481  in  the  system  with  unknown 
u  ,  V  ,  p  ,  and  p  . 


UI.  TECHNIQUE  OF  COMPUTATION 


The  solution  of  the  approximating  system  was  carried  out  numerically  with  a 
varying  niunber  of  intermediate  lines.  If  the  results  of  the  last  two  approxi¬ 
mations  coincided  with  sufficient  degree  of  accuracy,  this  was  taken  as  evidence 
of  the  practical  convergence  of  ths  computation.  As  already  noted,  integration 
of  the  system  was  carried  out  starting  at  the  axis  of  symmetry  $  •  Q  $  where 
Vfl  ■  -  0,  o>ir/2and  ^  andu|(0)(i  -  2, 3,...,  n)  are  the  N  unknown  parameters. 

The  principal  difficulties  in  the  calculations  lie  in  the  fact  that,  firstly,  it  is 
necessary  to  determine  these  N  unknown  parameters  from  the  conditions  on  the 
singular  line,  and  secondly  that  after  one  of  these  conditions  is  satisfied  at  each 
singular  point,  it  is  necessary  to  go  numerically  through  this  point  (in  the  solu¬ 
tion  it  is  necessary  to  go  through  all  the  N  singular  points  of  the  approximating 
system). 

At  each  singular  point,  at  least  one  of  the  equations  of  the  system  will  have  an 
indeterminate  right-hand  side  (0/0)  :  furthermore,  the  accuracy  of  the  calcvilation 
of  the  right-hand  sides  becomes  very  poor  in  the  vicinity  of  this  point,  and  there¬ 
fore  ordinary  methods  of  numerical  integration  become  unsuitable.  To  make  cal¬ 
culation  possible  there,  the  solution  of  the  system  in  the  vicinity  of  an  arbitrary 

nonsingular  (regular)  point(0  (0),  €(0) ,  e(0), . .  ,  ^  (0)  )  is  expanded  in  powers  of 

0  -  0(0) ,  and  in  the  vicinity  of  a  singular  point ( ,  «•  ,  a*  ,  . . . ,  ^)  it  is  expanded 
in  powers  of  0  -  0^(k  >  0,2,3,..;,  N).  The  convergence  of  such  series  was  de¬ 
monstrated  elsewhere  (15), 

In  the  calculation,  seversd  steps  ahead  of  the  singular  point  on  the  integral 
curve  (which  is  plotted  In  the  ordinary  manner)  a  regular  point  is  selected 
(0  m  0^^),  (  >  A  ■  0^^^} ,  **id  series  are  constructed  in  the  vicinity  of 

this  regular  point.  Tnese  series  are  used  to  extend  the  solution  to  th«^  singular 
point,  then  by  variation  of  the  parameter!  at  k  -  0  the  parameter  isf^  ,  and 
atk  >  2,3,..>,Nthe  parameter  is  s^0)lthe  required  condition 

2  2 

E||  ■  0  whes  -  (y  -  1  +  2s^)/(y'f  1) 

becomes  satisfied  at  the  k'th  singular  point.  Thus  we  determine  approximately 
the  value  of  the  parameter  and  the  value  of  the  quantities  ^  c* ,  a*, .  .  .  ,  ,  from 

which  we  then  build  up  the  series  in  the  vicinity  of  the  singular  point .  Using  the 
condition  that  atd  *  tfi^^these  series  must  Join  the  solution  obtained  in  the 
oridnary  manner,  we  refine  by  the  method  of  iteration  the  values  of  the  para¬ 
meter  and  the  values  of  the  quantities  <•,  o*,  .  .  .  ,  ;  the  number  of  unknows 

in  the  iterations  equals  the  number  of  conditions. 

In  the  calculation  of  the  flow  around  a  circular  cylinder,  both  types  of  series 
were  constructed  for  N  - 1,2,3  .  Examples  of  such  series  are  given  below  for 
the  first  approximation.  Furthermore,  it  is  also  known  that  if  the  singular 
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point  is  a  "saddle,"  then,  as  the  distance  from  the  saddle  increases,  all  the 
integral  curves  converge  to  two  integral  curves  that  pass  through  the  singular 
point  (in  our  problem  only  one  such  curve  is  joined  to  the  solution  ahead  of  the 
singular  point).  Using  these  series  and  the  property  of  convergence  of  the 
integral  curves,  we  succeeded  in  devising  technically  simple  ways  of  passing 
through  the  singular  points  of  the  approximating  system. 


The  system  of  ordinary  differential  equations  was  integrated  numerically  either 
by  the  Euler  method  with  iteration  or  by  the  Runge-Kutta  method.  The  inte¬ 
gration  was  in  steps  of  A9  =  2~7to  2*8;  the  functions  d  and  d  In  ^/d0  on  the 
wave  were  fixed  in  steps  of  Ad  ■  2-5  to  2~^,  The  values  of  the  functions 
d  In  ^ 

di  (^j)  and— -  were  found  on  the  intermediate  lines  by  quadratic  inter- 

d^* 

polation  from  the  values  of  the  corresponding  functions  on  the  wave.  The  series 
for  passing  through  the  singular  points  were  constructed  in  all  approximations. 
For  example,  for  =  5.  0  the  series  in  the  first  approximation ( N  -  1 )  were 
of  the  following  form: 


In  the  vicinity  of  the  regular  point  >  0.812)0) 


*  -  0.72151  ♦  0.74X2h  +  0.7)138h2  +  0.70674h5 
X  -  0.40377  +  0.44249h  -  0.3940)b^  +  0.)7399h5 
▼O  -0.39900  *  0.40036b  *  0.17)163b2  >  0.26982b) 


and  in  the  vicinity  of  the  singular  point  (  6g  >  0. 85791) 

<  -  0.7)699  ♦  0.81704b  0.7420)b2 

X  -  0.42264  0.48)29b  -  0.3)7278^ 

vg  -  0.4082)  *  0.37862b  0.1)18)08^ 


where 


X  »  fr/2  -  a  sod  h  ■  d  -  ^  . 


The  process  of  selecting  two  parameters  from  the  conditions  at  the  singular 
point  was  mechanized. 


IV.  RESULTS  OF  CALCULATIONS 


The  calculations  for  the  flow  past  a  circular  cylinder  were  carried  out  for 

«  2. 13,  2.5,  and  4.0  at  n  -  1  &nd  2,  and  for  M„  =  3.  0  and  5.  0  at  N  =1,  Z, 
and  3.  Some  of  the  results  of  these  calculations  are  presented  here  in  graphs 
and  in  tables. 

Figure  3  shows,  for  the  case  M„  =  3r  0  a^d  for  the  third  approximation  (N  =  3), 
the  shock  wave,  the  sonic  line  (  M  =  1),  and  the  limiting  characteristics  that 
bound  the  minimum  region  of  infltience.  The  region  of  influence  will  have  a 
similar  shape  .also  for  M„  -  Z,  5,  4.  0,  and  5.  0.  Figure  4  shows,  the  dependence 
of  (Q  on  found  in  various  approxir'iations.  Figures  5  2md  6  show  how  the 
distribution  of  the  pressure  varies  over  the  boCy,  and  also  the  shape  and  posi¬ 
tion  of  the  shock  wave  and  of  the  sonic  line  as  M,,  increases  from  2.  13  to  5.0 
[  here  pg(0)/pQ(O)  denotes  the  ratio  of  the  pressure  on  the  body  to  the  pressure 
at  the  stagnation  point  ]  .  The  results  are  shown  for  the  highest  approximation 
for  each  value  of  M„  .  For  M„  -  2. 13  and  2.  5,  no  sonic  lines  were  plotted  since 
the  calculations  for  these  Mach  niunbers  were  carried  out  only  for  N  =  1  and  2, 
and  the  accuracy  of  the  determination  of  the  sonic  lines  was  found  to  be  in- 
si^ficient. 

Figures  4, 7,8,  9^  and  10  illustrate  the  convergence  of  the  method  after  various 
approximations  for  different  values  of  .  Figure  4  shows  the  convergence  of 
the  method  in  finding  the  distance  ^  ,  while  figures  7  and  8  show  the  convergence 
in  finding  the  pressure  on  the  body  and  on  the  shock  wave  [p.  (d)/  pj  (0)  denotes  the 
ratio  of  the  pressure  on  the  wave  to^the  pressure  at  its  head]  .  Figures  9  and  10 
give  an  idea  of  the  accuracy  of  the  calculation  of  the  shape  and  position  of  the 
shock  wave  and  sonic  line  in  various  approximations  for  m„  -3.0  and  M,  =  5.  0. 

Most  papers,  both  here  and  abroad,  have  been  devoted  to  a  determination  of  the 
distance  from  the  surface  of  the  body  to  the  wav'.^  along  the  syiTunetry  axis. 
Figure  11  shows  a  comparison  of  the  results  obtained  experimentally  and  by 
calculations.  It  is  seen  that  our  results  are  in  good  agreement  with  the  experi¬ 
mental  data  (for  example,  for  M„  ■  4.0  Kim  (17)  gives  a  value  cq  s  0.  54,  while 
we  gat  <0  a  0.546). 

Figures  12,  13,  and  14  give,  for  various  values  of  M„  •  a  comparison  with  experi¬ 
ment  of  the  pressure  distribution  over  the  body.  The  experiments  for  M,.  s  2.  5 
and  3.  0  were  carried  out  by  the  scientific  assistant  G.  M.  Ryabinkov,  while  the 
experiments  for  a  4.  0  were  carried  out  by  Kim  (17). 

Figure  15  shows,  for  M„s  4.0,  the  shock  waves  obtained  by  the  method  of  integral 
relations  (N  a  2),  using  the  theory  of  K.  Tamada  (8)  and  the  experimentally  ob¬ 
tained  points  (17).  From  all  these  comparisons,  we  see  that  the  distribution  of 
pressure  over  the  body  and  the  shape  and  position  of  the  shock  wave  are  in  good 
agreement  with  experiment. 
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We  now  compare  the  sonic  lines.  Firstly,  Kim  (17)  in  his  experimental  work 
indicated  that  the  sonic  line  is  convex  at  M.,  =  4.0.  He  wrote:  "The  sonic  line 
is  bent  inside  near  the  wave.  It  is  not  clear  whether  this  is  a  very  important 
property  of  the  phenomenon,  or  whether  it  is  caused  by  the  fact  that  the  shock 
wave  is  somewhat  smeared  out,  a  fact  resulting  from  the  action  of  the  boundary 
layers  on  the  window  of  the  tube.  "  Thus,  the  result  we  obtained  (Fig.  6)  is  in 
qualitative  agreement  with  Kim's  result.  Furthermore,  in  (15)  we  obtained 
correct  expressions  for  the  angles  between  the  lines  w  =  const  and  the 
direction  of  the  incident  stream  in  the  plane  case.  It  was  shown  there  that  on 
the  shock  wave  these  angles  depend  only  on  the  magnitude  of  the  velocity  w  and 
on  the  Mach  number  M.,  of  the  incident  stream.  The  magnitudes  of  these  angles 
for  different  values  w  =  const  (  or  M  =  const)  and  M.,  are  shown  in  figure  17. 

For  the  sonic  line  (  M  =1),  we  have  by  calculaticn  the  following  values  for  these 
angles:  when  =  3.0,  we  have  V  20*  (  N  =  2)  and  =  1-2*  (N  =  3);  when 
=  4,'^B  -  21-22*  (  N  =  2),  and  when  =  5.0  we  get<l>B  S'  27-28*  (  N  =  2). 

The  exact  values  of  the  angle  for  the  same  values  of  are  3.  312*  ,  19.  840*  , 
and  26.  306*  ,  respectively.  Thus,  for  small  values  of  the  accuracy  with 
M  =  2  is  inadequate  for  plotting  sonic  lines,  and  when  M,.  s  4.  0  and  5.  0  even  the 
second  approximation  gives  sufficiently  accurate  results. 

A  Busemann  suggested  in  his  paper  (6)  that  the  sonic  point  on  the  body  coincides 
with  the  point  where  the  slope  of  the  surface  equals  the  critical  angle  of  wedge 
or  cone  (for  axially  symmetrical  flow),  at  which  the  frontal  shock  wave  becomes 
detached.  The  positions  of  the  sonic  points,  determined  from  this  assumption, 
are  compared  in  figure  16  with  our  own  results.  The  same  figure  shows  the 
experimental  points,  obtained  by  D.  Holder  and  A.  Chinneck  (18),  for  M„  =  1.42, 
1.60,  and  1.79.  It  is  seen  from  the  graph  that  the  sonic  points,  found  by  the 
method  of  A,  Busemann,  are  located  at  much  higher  points  on  the  body  than  the 
points  obtained  either  by  the  method  of  integral  relations  or  experimentally. 
Furthermore,  it  follows  from  the  graph  that  as  the  Mach  number  increases 
the  sonic  point  shifts  on  the  body  towards  the  stagnation  point,  a  fact  to  which 
Holder  and  Chinneck  called  attention. 

An  analysis  of  the  calculation  and  a  comparison  of  the  results  obtained  with 
experiment  show  that  in  the  vicinity  of  the  axis  of  symmetry  even  the  first 
approximation  gives  fair  results.  The  distribution  of  pressure  on  the  body  and 
on  the  wave  and  the  shape  and  position  of  the  shock  wave  are  all  determined  with 
sufficient  accuracy  over  the  entire  region  of  integration  by  means  of  the  second 
approximation;  to  construct  the  ionic  lines  at  <  3.  0  to  3.  5  we  must  have  at 
least  three  approximations,  while  for  M.  >  3.  5  two  appro.ximations  are  sufficient. 

The  computed  data  make  it  possibb  to  determine  the  distribution  of  the  pressure 
on  the  surface  and  outside  the  body,  the  sonic  line,  the  characteristics,  the 
position  and  shape  of  the  shock  wave,  and  so  on. 

In  conclusion,  we  tabulate  the  results  of  the  calculations. 
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In  Table  I,  we  give  the  values  of  the  following  quantities  as  functions  of  the 
angle  6  : 

(  --  distance  from  the  body  to  the  shock  wave  along  $  =  const 

X  m/2  —  o  inclination  of  the  wave  to  the  vertical 

entropy  fvinction  (the  entropy  is  S«(y-i)  c,ln  ) 

d  liD 

- i -  ~  derivative  of  the  entropy  function  with  respect  to  ^ 

Po(^)/Po(0)  ~~  ratio  of  the  pressure  on  the  body  to  the  pressure  at  the 

stagnation  point. 


Tables  II,  HI.  and  IV  list,  for  various  u  ,  the  values  of  the  flow  variables 

«• 

behind  the  shock  wave  at  five  points  on  the  ray  6  -  const  (the  distance  from  the 
body  to  the  wave  along  B  -  const*  is  divided  into  four  parts:) 


t  -  1 


s  0  (body):  0.25;. 0,50,  0.75,  and  1.00  (shock  wave) 


>•  velocity  components  along  r  and  9 


X,  p,  p  ><*  temperature,  density,  and  pressure 


^  —  stream  function 

S9  »3.5«  . 


Tables  V  and  VI  give  the  coordinates  of  the  sonic  line  and  of  the  limiting 
characteristic  for  various  values  of  . 


Table  VU  gives,  for  various  values  of  •  the  values  of  the  angles  9,^  and 

at  which  the  sonic  line  and  the  limiting  characteristic  intersect  the 
surface  of  the  body. 

In  all  the  tables,  the  angles  are  in  radians,  and  the  remaining  quantities  are  in 
dimensionless  form:  linear  dimensions  are  referred  to  the  radius  of  the  cylinder; 
the  velocity  is  referred  to  the  maximum  velocity:  the  density,  pressure,  and 
temperature  are  referred  to  the  corresponding  stagnation  values  ahead  of  the 
shock  wave. 
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0312 

0194 

395 

843  . 

23* 

196 

0533 

0625 

0288 

403 

830 

231 

193 

0550 

0938 

0385 

411 

829 

228 

189 

0^ 

1,1250 

0,0485 

0,418 

0,822 

0,224 

0,184 

0,0584 

1562 

0588 

425 

815 

221 

180 

OWl 

l«0.T5 

0,0000 

-0,168 

0,0000 

0,972 

0,805 

0,297 

0,0000 

0625 

—  165 

0295 

.  071 

805 

297 

00485 

1250 

—  162 

0501 

069 

805 

295 

00070 

1875 

—  150 

0887 

966 

303 

203 

0145 

2500 

—  154 

118 

962 

801 

200 

0104 

0,8125 

-0,140 

0,148 

0,956 

0,290 

0,286 

0,Q^?2 

8750 

—  141 

175 

949 

296  . 

281 

0291 

4375 

—  133 

202 

041 

293 

270 

0340 

5000 

—  123 

220 

932 

280 

270 

0380 

5625 

—  Ill 

255 

922 

285 

263 

0430 

•0,6250 

-0,0090 

0,280 

0,011 

0,281 

0,256 

d,0480 

6875 

-0850 

304 

890 

276 

248 

0530 

33 
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'1 


5-1,00 


1  -0,207 

0,0000 

0,056 

0.294 

0,281 

-  207 

0387 

056 

293 

280 

■  —  204 

0671 

953 

293 

279 

-  199 

0098 

950 

292 

278 

-  193 

131 

945 

202 

276 

-0,185 

0,162 

0,838 

0,291 

0,273 

—  175 

182 

881 

289 

269 

-  164 

222 

028 

288 

266 

-  151 

281 

014 

286 

261 

-  186 

278 

908 

284 

257 

-0,120 

0,808 

0,882 

0,282 

0,252 

-  M2 

828 

880 

280 

246 

-  0822 

842 

974 

270 

248 

-  0121 

188 

m 

277 

241 

-  0720 

184 

801 

276 

238 

-0,0818 

0,878 

0,888 

0,275 

0,235 

—  0888 

888 

860 

274  ■ 

233 

-  OMO 

884 

848 

272 

229 

-  0288 

404 

888 

271 

226 

-  0187 

418 

829 

260 

223 

-0,00408 

0,421 

0.622 

0,268 

0,220 

00708 

428 

818 

266 

217 

0208 

487 

806 

264 

214 

0881 

444 

801 

263 

210 

0463 

451 

794 

261 

207 

0,0688 

0,487 

0,787 

0,260 

0,204 

0732 

402 

780 

258 

201 

0,0000 

00650 

0130 

0195 

0261 

0,0327 

0393 

0461 

0529 

0598 

0,0668 

0739 

0778 

0815 

0652 

0,0890 

0919 

0965 

100 

104 

0,108 

112 

117 

121 

125 


34 


TABLE  m 


• 

M 

V 

T 

t 

p 

0,0000 

0,0000 

0,0000 

5-0,00 

1,000 

0,138 

0,138 

0,0000 

0025 

0000 

0282 

0,888 

138 

138 

0000 

1250 

0000 

0598 

996 

137 

137 

0000 

1875 

0000 

0823 

991 

135 

134 

0000 

2500 

0000 

125 

084 

133 

131 

0000 

0,3125 

0,0000 

0,159 

0,874 

0,130 

0,126 

0,0000 

3750 

0000 

192 

063 

126 

121 

0000 

4375 

0000 

224 

949 

121 

115 

0000 

5000 

0000 

256 

934 

117 

109 

0000 

5625 

0000 

287 

917 

111 

102 

0000 

0,6250 

0,0000 

0,317 

0,890 

0,106 

0,0956 

0,0000 

6875 

0000 

347 

879 

loa 

0883 

0000 

7500 

0000 

377 

857 

om 

0811 

0000 

8008 

0000 

400 

839 

0895 

0752 

0000 

8438 

(UXX) 

419 

824 

0855 

•  0704 

0000 

0,8750 

0,0000 

0,433 

0,811 

0,0823 

0,0668 

0,0000 

9375 

0000 

461 

787 

0763 

0601 

0000 

1,0000 

0000 

487 

762 

0704 

0536 

oooo 

0625 

0000 

512 

736 

0646 

0476 

0000 

1250 

0000 

537 

711 

0582 

0421 

oooo 

1,1875 

o;oooo 

0,561 

0,685 

0,0539 

0,0369 

0,0000 

0,0000 

-0,0564 

0,0000 

5-0,25 

0,808 

0,197 

0,136 

0,0000 

0625 

—  0563 

0907 

905 

136 

136 

000571 

1250 

^  0550 

0611 

003 

135 

135 

00114 

1875 

-  0548 

0012 

868 

184 

133 

00171 

2500 

-  0534 

120 

982 

132 

130 

00228 

0,3125 

-0,0513 

0,150 

o,r4 

0,130 

0,128 

0,00284 

8750 

-  04r 

170 

865 

ia 

122 

00338 

4375 

-  0455 

207 

954 

128 

118 

00382 

5000 

-  0417 

885 

642 

120 

113 

00445 

5625 

-  0672 

288 

828 

116 

107 

00486 

0,6250 

-0,0918 

0,218 

0,915 

0,111 

0,102 

*  0,00545 

9675 

-  0980 

815 

888 

107 

0066 

00593 

7500 

-  0188 

m 

868 

102 

0006 

00638 

•ns 

-  0181 

880 

660 

0868 

0860 

00674 

8438 

—  00754 

878 

806 

0857 

0821 

00703 

0,6750 

-0,00821 

0,886 

0,648 

o;oo32 

0,0702 

0,00728 

0375 

00608 

410 

881 

0885 

0730 

00762 

1,0000 

0162 

481 

818 

0697 

0681 

00790 

0625 

0272 

450 

795 

0790 

0628 

00884 

1250 

0881 

408 

778 

0744 

0579 

00666 

1,1875 

0,0520 

0,467 

0,780 

0,0699 

0,0531 

0,00606 
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TABLE  m  (Cont'd) 


€  -0.50 

0,0000  —0,106  0,0000  0.B88  0,134  0,133  0,000 

0625  —  106  0310  987  134  132  00115 

1250  —  135  0631  984  133  131  002^ 

1875  —  103  0931  980  132  130  00345 

2500  —  100  122  975  131  127  00460 

0,3125  —0,0963  0,150  0,968  0,129  0,125  0,00574 

3750  —  0911  178  959  127  122  00687 

4375  —  0848  206.  950  125  118  b07M 

5000  —  0773  233  939  122  115  00908 

5625  ^  0687  259  928  119  111  0101 

0,6250  —0,0588  0,284  0,915  0,117  0,107  0,0112 

6875  —  0477  309  902  113  102  0122 

7500  —  0352  332  888  110  0985  0133 

8008  —  0241  350  876  108  0949  0141 

8438  -  0140  364  *866  106  0919  0148 

0,8750  —0,00626  0,374  0,859  0,104  0,0897  0,0153 

9375  0,0102  394  844  101  0854  0162, 

1,0000  0281  412  829  0978  0811  0172 

0625  0474  428  814  0945  0769  0181 

1250  0680  443  798  0913  0729  0190 

1,1875  0,0900  0,457  0,7K  0,0879  0,0688  0,0199 

«-0,75 

0,00(Ki  I  -0,151  0.0000  0,977  0,131  |  0,128  0,0000 

0625  !  -  151  0830  VK  131  127  00173 

1250  -  149  0857  973  130  127  00348 

1875  -  146  0979  969  129  125  00622 

2500  -  141  129  983  129  124  00898 

0,3125  -0.184  0,180  0,958  0,128  0,122  0,00870 

3750  -  127  190  947  126  120  0104 

4378  —  117  220  987  12^  117  0121 

5000  —  107  .  248  928  124  115  0139 

8825  -  0940  270  914  -  122  112  0156 

0,6280  -0,0907  0,802  0,902  0.120  0,109  0,0173 

9878  -  0081  827  898  119  106  0190 

7800  —  0491  881  874  117  102  0207 

8008  -  0880  888  862  116  100  0221 

8438  —  0104  898  852  114  0977  0233 

0,8780  -0,00918  0,308  0,844  0,113  0,0961  0,02i2 

9378  0,0125  412  829  ill  0028  0259 

1,0000  0358  430  313  109  0895  0276 

0625  0006  445  797  107  0861  0294 

1280  0867  489  781  106  0828  0312 

1,1875  0,114  0,471  |  0,764  0,103  0,0794  |  0,0630 
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TABLE  m  (Cont'd) 


• 

II' 

• 

^  1 
1 

t 

p 

♦ 

0,0000 

-0,190 

0,0000 

6-1,00 

0,963 

0,126 

0,121 

0,0000 

0625 

—  180 

0339 

962 

126 

121 

00233 

1250 

—  187 

0690 

960 

126 

121 

00467 

1875 

—  182 

105 

955 

126 

120 

00701 

2500 

—  175 

143 

948 

125 

119 

00037 

0,3125 

—0,167 

0,179 

0,939 

0,125 

0,117 

0,0117 

3750 

—  156 

215 

928 

124 

116 

0141 

4375 

—  lU 

249 

916 

124 

114 

0164 

5000 

—  130 

282 

903 

123 

111 

0189 

5625 

—  114 

313 

888 

123 

100 

0213 

0,6250 

-0,0976 

0,343 

0,872 

0,122 

0,106 

0,0237 

6875 

—  0783 

371 

855 

121 

104 

0262 

7500 

—  0577 

307 

838 

120 

101 

0287 

8008 

—  0307 

417 

823 

119 

0087 

8438 

—  0238 

434 

810 

119 

0066 

0326 

0,8750 

—0,0118 

0,445 

0,801 

0,118 

0,0050 

0,0339 

9375 

0129 

466 

782 

117 

0018 

0386 

1,0000 

0302 

485 

762 

116 

0686 

0303 

0625 

0606 

502 

743 

114 

0854 

0422 

1250 

0050 

517 

723 

113 

0821 

0451 

1,1875 

0,124 

0,529 

0,704 

0,112 

0,0789 

0,0481 

TABLE  IV 

• 

< 

9 

r 

f 

1  ' 

♦ 

1-0,00 

0,0000 

0,0000 

0,0000 

1,000 

0,«17 

0,M17 

0,0000 

0818 

0800 

OHO 

0,9H 

M15 

0615 

oooo 

un 

omo 

0817 

OH 

Mil 

0</09 

oooo 

wn 

0000 

oon 

OH 

OOM 

J607 

oooo 

MO 

oooo 

128 

OH 

COM 

OMt 

oooo 

0,8128 

0,0000 

0,1H 

0,078 

O.M70 

0,OM1 

0,0000 

mm 

0000 

IH 

081 

OOM 

0537 

oooo 

4878 

oooo 

2H 

047 

06M 

0510 

oooo 

8000 

oooo 

HI 

HI 

M17 

0481 

oooo 

mm 

oooo 

2R 

014 

OHO 

0451 

oooo 

0,0280 

0,0000 

0,8» 

O,0H 

0,04M 

C,0410 

0,0000 

0078 

oooo 

H4 

074 

0441 

MM 

oooo 

7800 

oooo 

m 

OH 

0414 

OSH 

oooo 

0128 

oooo 

411 

OH 

0388 

0322 

oooo 

mm 

oooo 

441 

006 

OOM 

02M 

,A)00 

37 


TABLX  IV  (Cont'd) 


• 

m 

-  t 

r 

t 

P- 

v 

♦ 

0,0375 

0,0000 

0,469 

0.779 

0.0331 

9,0258 

0.0090 

1,0000 

0060 

496 

754 

0SD4 

022» 

0090 

0625 

0000 

522 

121 

0278 

0202 

oo« 

1250 

0000 

547 

700 

0253 

0177 

0099 

1875  : 

0000 

570, 

674 

0230 

0155  1 

0090 

1-0,25 


6,0000 

—0,0510 

0,0000 

0,997 

0,0611 

0,0609 

0,0000 

0625 

—  0520 

0316 

996 

0610 

0607 

OP0230 

1290 

—  0514 

9030 

993 

0605 

0601 

000460 

1875 

—  0504 

0039 

988 

0598 

0592 

000668 

2500 

—  0480 

124 

962 

0589 

0578 

000915 

0,8125 

—0.0460 

0,154 

0,973 

0,0577 

0,0562 

0,00113 

3750 

—  0444 

184 

964 

0563 

0543 

00135 

4376 

—  0413 

213 

952 

0547 

0521 

00156 

5000 

—  0376 

242 

939 

0530 

0498 

00177 

5625 

—  0333 

270 

925 

0511 

•  0473 

00197 

0,6250 

-0,02tt 

0,208 

0,9*0 

0,0491 

0,0447 

0.00216 

68TS 

—  0224 

325 

803 

0470 

0420 

00235 

7110 

—  0158 

359 

876 

C448 

0393 

00252 

•m 

—  005«2 

374 

839 

0426 

0366 

00269 

1710 

—  000141 

309 

910 

0408 

0339 

00284 

0,9375 

O.OOI59 

0,421 

0,122 

6,0381 

9,0313 

0,00298 

1,0000 

0190 

443 

m 

0868 

0288 

00311 

005 

0869 

455 

m 

0336 

6264 

00323 

1250 

3417 

48it 

705 

0315 

0241 

06334 

1875 

05a 

500 

747 

0294 

0220 

00344' 

9,0000 

-0.0097 

0.0000 

0125 

-  ilN 

0112 

1250 

-  om 

0055 

1075 

-  0002 

0105 

1500 

-  <P0 

181 

0,5125 

-o.oov 

0,150 

9750 

-  0085 

100 

4575 

-  0778 

214 

SOUO 

-  0702 

243 

5625  J 

—  0517 

270 

0,0250 

-0.0619 

0,296 

0075 

—  0408 

322 

7500 

—  0285 

1  SM 

8125 

—  0147 

368 

.  5750 

o,oou 

390 

0.990 

0,0002 

0,0596 

909 

0000 

0594  : 

106 

0597 

0588 

991 

0591 

0580 

975 

0584 

0570 

0.907 

0,0576 

0,0558 

059 

0567 

0543 

947 

0555 

0526 

936 

0543 

0509 

923 

o^ao 

0489  ' 

0,909 

0,0516 

0,0469  1 

894 

0502 

0449  : 

879 

0487 

0428  ! 

863 

0472 

0407 

847 

0456 

1 

0386  ' 

0,0000 

000460 

000935 

00139 

00186 

0,00232 

00277 

00321 

00365 

00408 

0,00450 

00402 

00532 

00571 

00600 


38. 


TABLE  IV  (Cont'd) 


» 

m 

9 

T 

r 

p 

♦ 

0,S375 

0,0168 

0,410 

0.831 

0,0441 

0,0366 

0,00646 

1.0000 

0347 

428 

815 

0425 

0346 

00682 

0825 

0540 

445 

798 

0409 

0327 

00716 

1250 

0747 

461 

781 

0394 

0308 

00750 

1875 

0068 

474 

765 

0378 

0289 

00782 

£-0,75 

0.0000 

-0,143 

0,0000 

0,979 

0,0587 

0,0575 

0.0000 

0625 

—  142 

0348 

978 

0586 

0573 

000709 

1250 

—  140 

0603 

975 

0584 

0569 

00141 

1875 

—  137 

103 

970 

0581 

0564 

00212 

2500 

-  132 

136 

963 

0577 

0356 

00283 

0.3125 

-0,125 

0,169 

0,955 

0,0572 

0,0547 

0,00354 

3750 

—  118 

201 

945 

0567 

0536 

00424 

4375 

—  108 

232 

934 

0561 

0524 

00494 

5000 

—  oon 

262 

921 

0555 

0511 

00564 

5625 

—  0852 

291 

907 

0547 

0497 

00633 

0.6250  ' 

-0,0711 

0,319 

0,893 

0,0540 

0,0482 

0,00702 

6875 

-  0554 

345 

877 

0532 

0467 

00771 

75<10 

-  0381 

370 

0451 

8125 

-  0191 

583 

844 

0515 

0435 

8750 

0,00144 

414 

828 

0007 

0420 

00975 

0.9375 

0,0235 

0,434 

0,810 

0,0488 

0,0404 

0,0104 

1,0000 

0472 

452 

785 

0480 

0318 

0111 

0625 

0725 

408 

775 

0481 

0372 

0117 

1250 

0090 

412 

757 

0471 

0857 

0124 

1875 

127 

488 

758 

0482 

0341 

0181 

1-1,00 

0,0000 

-0,182 

o.om 

0,960 

0,0587 

0,0548 

0,0000 

0825  i 

-  IM 

0884 

015 

0686 

0547 

000958 

1250  ' 

-  178 

0741 

912 

0566 

0645 

00191 

1875 

-  178, 

118 

980 

0665 

0641 

00286 

25l« 

-  188 

158 

848 

0664 

0536 

00384 

0,3125 

-0.158' 

0,102 

0,887 

0,0658 

0,0628 

0,00481 

3750  > 

-  147 

230 

828 

0582 

0520 

00578 

4375 

-  134 

200 

810 

0680 

0610 

00675 

9INI0 

-  120 

300 

885 

0558 

0409 

00773 

5625 

—  108' 

)33 

877 

0555 

0487 

00672 

0,6250 

-0,0058 

0,364 

0,850 

0,0553 

0,0475 

0,00071 

6875 

-  0850 

394 

840 

0550 

0462 

0107 

7.500 

-  0445 

422 

810 

0546 

0448 

0117 

8125 

-  0213 

448 

798 

0543 

0433 

0127 

8750 

0,00317 

472 

776 

0530 

0419 

0138 
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TABJLE  IV  (Cont'd) 


• 

1 

1  * 

1 

9' 

r 

r 

P 

♦ 

0,8375 

0,0291 

0,494 

0,754 

0,0535 

0,0404 

0,0148 

1,0000 

0563 

514 

731 

0531 

0388 

0159 

0625 

0849 

532 

709  . 

0526 

0873 

0170 

1250 

114 

548 

686 

0521 

0357 

0182 

1675 

145 

561 

664 

0516 

0342 

0194 

TABLE  V 


- 

3,f 

-«.o 

-&.0 

y 

•—  X 

y 

—JT 

+y 

0,664 

0,748 

0,681 

0,732 

0,693 

0,721 

666 

708 

082 

no 

710 

798 

668 

852 

704 

843 

724 

867 

600 

911 

716 

014 

745 

951 

670 

974 

731 

995 

786 

1,055 

0,671 

1,045 

0,757 

1,102 

0,815* 

1,097 

on 

135 

ooo* 

217 

830 

113 

082 

240 

006 

280 

900 

160 

720 

840 

005 

314 

944 

167 

m* 

443 

1,020 

314 

976 

168 

0,841 

1,510 

1,000 

1,302 

1,006 

1,167 

840 

800 

148 

204 

084 

162 

1.080 

008 

108 

270 

007 

148 

117 

020 

101 

200 

188 

128 

270 

000 

810 

222 

247 

068 

TABLE  VI 


M 

Mg,  -  e,4 

— JT 

y 

-  1  7 

-*  1  y 

! 


0,575 

0,818,  V 

0,628 

0,778  1 

0,649 

0,761 

576 

834 

643 

837 

656 

784 

584 

904 

661 

901 

667 

816 

503 

972 

182 

964 

678 

947 

605 

1,030 

704 

1,024 

690 

on 

40 
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TABLE  VI  (Cont'd) 


-  3.0  ; 

it(«.  -  4.0  1 

&.0 

1 

1 

y  1 

-  1 

y 

— Jf 

y 

0,620 

1,101 

0,717 

1,054 

0,702 

0,907 

638 

163 

730 

082 

716 

937 

660 

222 

744 
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